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Abstract 

We study the time and temperature dependent correlation functions for an impenetrable bose gas 
with open boundary conditions. We derive the Fredholm determinant formulae for the correlation 
functions, by means of the Bethe Ansatz. In the case of time independent ground state, our Fredholm 
determinant formulae degenerate to the one which have been obtained by the help of fermions 



1 Introduction 



In the standard treatment of quantum integrable models, one starts with a finite box and impose 
periodic boundary conditions, in order to ensure integrabihty. Recently, there has been increasing 
interest in exploring other possible boundary conditions compatible with integrabihty. These other 
possible boundary conditions are called "open boundary conditions" . 

With open boundary conditions, the works on the two dimensional Ising model are among the 
earliest. By the help of graph theoretical approach, B.M. McCoy and T.T. Wu Q studied the two 
dimensional Ising model with open boundary conditions. They calculated the local magnetizations. 
E.K. Sklyanin Q began the Bethe Ansatz approach to open boundary problems. M. Jimbo et.al. 
10] calculated correlation functions of local operators for the antiferromagnetic XXZ chains with 
open boundary conditions, using Sklyanin's algebraic Bethe Ansatz framework and representation 
theoretical approach invented by Kyoto school [HQ. T. Kojima H] studied the ground state 
correlation functions for an impenetrable Bose gas with open boundary conditions : 

(V^(xi)Vt(x2)). (1) 

He derived the Fredholm minor determinant representations for the ground state correlation func- 
tions by the help of fermions, which have the integral kernel : 

sin(A - //) ^ sin(A + n) 
X- H X + n ' 

The integral intervals depend on the space parameter xi,X2- In this paper we study an impenetrable 
Bose gas with open boundary conditions. We are interested in the finite-temperature dynamical 
correlation functions : 

{i>ixuO)iP\x2,t))T. (3) 

We derive the Fredholm determinant representations for the dynamical correlation functions by the 
coordinate Bethe Ansatz, which have the integral kernel : 

L{X,fi) ^ L{X,-fj.) 



X — n X + ij, ^ 



where we have used 



L(A,/i) = e**^'sin(xi(A-/i)) + e**^%in(x2(A-^)) (5) 

2 f 1 1 \ • 2 

H — / e**'' sin((s - sin((s - X)x2)ds. 

vr J-oo \s — u s — XJ 
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The mesure of integral depends on the temperature parameter T > 0. For the special case that 
time t = and temperature T = 0, our formulas degenerate to the one which has been obtained 
||^]. In the case of periodic bounadary conditions, V.E. Korepin and N.A. Slavnov calculated the 
dynamical correlation functions for an impenetrable Bose gas by the corrdinate Bethe Ansatz. 

Now a few words about the organization of the paper. In section 2 we formulate the problem. 
In section 3 we obtain the determinant formulae for the field form factors. In section 4 we obtain 
the Fredholm determinant representation for the dynamical correlation functions. In section 5 we 
consider the special case that time t = and derive the Fredholm minor determinant representa- 
tions for the finite-temperature fields correlation functions. We show that our Fredholm formulae 
degenerate to the one which has been obtained |||]. 

2 Formulation 

The purpose of this section is to formulate the problem. The Hamiltonian of our model is given by 

H = £ dx (d^i^^d^ijj + cip^^^i)^ - hi)^^^ + ho (^t (0)^(0) - ^/;t(L)V'(L)) . (6) 

Here the fields ip{x) and ip'^{x) {x € R) are canonical Bose fields given by 

mx),i^Hy)]=6{x-y), [^(x), = [V'^(x), ^t(y)] = q, (x,yeR), (7) 

and L > is the size of box. The parameters h > and Hq represent the chemical potential and 
the boundary chemical potential respectively. We only consider the case of the coupling constant 
c = oo, so-called "impenetrable case". The Hamiltonian H acts on the Fock space of the Bose fields 
defined by the following relations between the Fock vacuum |0) and the Bose fields : 

(0|Vt(x) = 0, ^(x)|0) = 0, (0|0) = 1. (8) 

A A^-particle state vector I^'a?) is given by 

|*iv) = / dzi--- [ dzNTpNizi,--- ,ZN)i^Hzi)---i^HzN)\0), (9) 

Jo Jo 

where the integrand 'ip]\f{zi, ■ ■ ■ , zn) is a C-valued function. The eigenvector problem : H\^]\f) = 
EnI'^n), (En G R), is equivalent to the quantum mechanics problem defined by the following four 
conditions of the integrand function ipN^zi, • • • , zn)- 
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1. The wave function V^at = tpN{zi, • " ^ ^n) satisfies tlie free-particle Schrodinger equation in the 
case of variables < Zi ^ Zj < L : 

^ f d Y 

~Ylyg^.j i^N{zi,- ■ ■ ,zn) = En ■ i1^n{zi,- ■ ■ ,zn), {0 < Zi Zj < L, En e K). (10) 

2. The wave function -i/'Af is symmetric with respect to the variables : 

ifjNizi,--- ,Zn) ='>pN{Za{l),--- ,Za(N)), {cT ^ Sn)- (H) 

3. The wave function tpN satisfies the integrable open boundary conditions : 

= 0, U = !,■■■, N). (12) 



TT hollpN 



9zj 



0, f^ + ^ojV'JV 



4. The wave function ipN vanishes whenever the corrdinates coincide : 

IpNizi, ■ ■ ■ ,Zi,- ■ ■ ,Zj, - ■ ■ , ^Ar)!,.^,- = 0- (13) 

The wave functions '07V which satisfy the above four conditions were constructed [|l|. They are 
parametrized by the spectral parameters. 

iPn{zi, - ■ ■ ,zn\\i, - ■ ■ ,^n) 
= Cons. sgn{zj — Zk) det {XjCos{XjZk) + hQsm{XjZk))- (14) 

l<j<k<N ^<j,k<N 

Here the function sgn(x) = -j^ and the spectral parameters < Ai < A2 < • • • < Xn are determined 
by the so-called Bethe Ansatz equations : 

^i = pj' (^j-eN, J = 1,2, •••,]¥). (15) 

The constant factor "Cons." is deteminend by 

(^'^(Ai, • • • , AAr)|^'^(Ai, • • • , Xn)) = {2Lf. (16) 

The eigenvalue En{{X}) : 

H\^n{Xi, • • • , XN))e = En{{X})\^n{Xi, XN))e, (17) 

is given by 

TV 

ENi{X}) = T.i>^]-h). (18) 
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We assume that the set {|^'Ar(Ai, • • • , AAr))}aii{A}jv a^gn is a basis of physical space of this model. 
Here the index all{A}Ar represents all the solutions of the Bethe Ansatz equations (|l5|). This type 
assumption is usually called "Bethe Ansatz" . The following lemma is a foundation of our analysis. 

Lemma 2.1 // the boundary condition /iq takes the special value ho = 0,oo, the eigenvectors 
|^'({A})) satisfy orthogonality relations. 

N 

(*^(Ai,---,Aw)|^'^(^i,---,^^)) = (2L)^n5A„M.' (/io = 0,cx)). (19) 
Here 5x^^ is Kronecker Delta. 

To prove the above lemma, we have used the Bethe- Ansatz equations of the spectral parameters. 
In the sequel we use the orthogonality relations of the eigenstates, therefore we concentrate our 
attentions to the case of the special bounadary conditions : Hq = 0, oo. The boundary conditions 
/iQ = and ho = oo are called Neumann, Dirichlet, respectively. In the sequel we use the following 
abberivations. 

|^'Ar(Ai, • • • , Xn))+ for Neumann, |^'Ar(Ai, • • • , Xn))~ for Dirichlet. (20) 
The constant ^^Cons.^^ is given by 

(rijLi Aj^ , for Neumann, 



Cons. = < 



2N 

Vm(S)''' for Dirichlet. 



From the orthogonal relations (^) and the so-called "Bethe Ansatz" , we arrive at the completeness 
relations. 

The Bose fields 'ip{x,t),'il)^ {x,t) are developed by the time t by 

idtiJ = [^,H], idt^^ = [ijlH]. (23) 
More explicitly the time dependence of the Bose fields are written by 

iPix, t) = e'^'i^{x)e-'^\ V^(x, t) = e'^'iP\x)e-'^K (24) 

In this paper we are interested in the dynamical correlation functions {4^{xi,ti)'ilj^ {x2,t2))e,T de- 
fined by the following way. For the nonzero temperature T > 0, the dynamical correlation functions 
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for N state are defined by the summation of the every states : 

{Hx^,t^)^|;Hx2,t2))e,N,T = { E exp (-:^^^^) 

^ / giv({A}) \ e(^jv(Ai,---,Ajv)|V;(xi,ti)V.t(a:2,t2)|^jv(Ai,---,Ajv)), | 
.aiiTTu''''^^ r J ,(*^(Ai,---,A;v)|*iv(Ai,---,Aiv))e J'^ ^ 

where the index e = it represents the boundary conditions. For the ground state case T = 0, the 
dynamical correlation functions for AT state are defined by the vacuume expectation value of the 
ground state : 

. ^ lU + ^\ e(*Af(Ai,---,Aiv)|V'(xi,ii)V'Ha;2,i2)|*iv(Ai,---,Aiv))e 

where the spectral parameters (Ai, • • • ,Xn) are given by 

A - = / 5^-^ ~ Neumann, ^^^^ 

I J J, for Dirichlet. 

In this paper we are interested in the thermodynamics of the correlation functions. For the nonzero 
temperature T > 0, the dynamical correlation functions in the thermodynamic limit are defined by 



{ll^{xi,ti)'ll^^{X2,t2))e,T = Jim {il;{xi,ti)'il;^X2,t2))e,N,T- (28) 

Here the density D{T) = is given by 

D{T) = 3 / ■&{X)dX, (29) 
where the Fermi weight ■i?(A) is given by 



TT JO 



For the ground state T = 0, the dynamical correlation function in the thermodynamic limit is defined 

by 

{lp{xi,ti)'tp^{x2,t2))efl= {Tp{xi,ti)lJj^X2,t2))e,N,0- (31) 

^=0(0) 

Here the density D{0) = can be chosen arbitrary. In this paper wc give the Fredholm determinant 
representations for the dynamical correlation functions {il){xi,ti)ip'' {x2,t2))e,T, {T > 0,e = ±). 
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3 Form Factors 

The purpose of this section is to derive the determinant formulae for the form factors. First we 
prepare a lemma. 

Lemma 3.1 For the sequences {fj^k}j=i,---,N+i, k=i,---.N o,nd {gij}j=i^...^Ar_|_i, the following holds : 



5^ sgna/^(^+i) n 9o(3),3 = [fN+i + -Q^j ^ fe.fc - afj ■ gN+i,k) 



(32) 

0=0 



Proof. Consider the coset decomposition : 

Sn+1 = Sn{N + 1) U Sn{N) • (A^, iV + 1) U • • • U Sn{1) ■il,N + 1), 
where S^ij) is permutations of (1, • • • , j — 1, + j + 1, • • • , A^). Rewrite the left side of the equation 



(3.1) with respect to the coset decomposition : 



N N 

(L.H.S.) = J29j H '^Sn(r • (j, iV + 1)) [| f^(^k),k ■ fr(N+i),j 

N 

+ QN+i sgn(r-(iV + l,iV + l))n/r(fc),fc = (^-^^.5-) (33) 

Te5jv(Af+l) k=l 



Now let us consider the field form factor : 



/o JO 
1 



(1 + '5ai,o)(1 + ^fiifi) ct&Sn+i 



W dz sgn(z - x)(e-*^-(^)^ + ee^^-0)^)(e^^^^ + ee-'^^J^)! . 



To derive the third line, we have used a simple fact : 



Using lemma 3.1, we arrive at the determinant formulae for the form factors. 



□ 



e(*Ar+l(Ai, • • • , \N+l)\lp\x)\^!N{lJil, • • • , /UAf))e 

y/N + 1 / dzi--- dzNtp*N+i{zi,- ■ ■ ,ZN,x\\i,- ■ ■ ,Xn+i) (34) 



sgna(e-^^-(^+i)^ + ee^^-(^+i)^) (35) 



N N 

J2 sgnarn/a(i),r(i) = A^! I] sgnan/a(i),i- (36) 

o-,Te5]v j=i o-eSjv j=i 



Lemma 3.2 The field form factors have the determinant formula. 

e(^Ar+l(Ai, • • • , AAr+i)|V'"l'(x)|^'Ar(/ii, • • ■ ,fJ.N))e 

= [C^{x\Xn+i) + tt) det {I^{x\Xj,Hk) - aCe{x\Xj)I^{x\XN+i,fJ.k)) ■ (37) 
V oaj i<j,k<N 

Here we have ued 

a(x|A) = —L={e-'^- + ee'^-), (38) 

1 f 4 

/,(x|A,/i) = , { -. sin(x(A-;u)) (39) 

V(i + M(i + W ^^"^ 

4 1 
+ e-r— — sm(3;(A + /i)) - 2L{5x^^ + eSxfld^fl) \ . 

A + fl J 



We can write the field form factors without using integrals. From the relation ^/^^ {x, t) = e*^*'0^ {x)^ *^*) 
the dynamical form factors are given by 

e(^'Ar+i(Ai, • • • , AAr+i)|V'"'"(a:,t)|^'Ar(/Ui, • • • ,/UAr)), 
= expjii A^ -^^^ij I ,(*iv+i(Ai,---,Aiv+i)|V'"^(x)|^'iv(Mi,---,m)).. (40) 

4 Correlation Functions 

The purpose of this section is to derive the Fredholm determinant formulas of the dynamical cor- 
relation functions {ip{xi,ti)ip'^ {x2,t2))e,T- First we consider the vacuum expectation values of fields 



operators. Using the completeness relation (|22D , the vacuum expectation values of two fields are 
given by 

g(^'iv(/^l, • • ■ ,IJ-N)\lPixi,ti)lp^X2,t2)\'^NifJ'l, ■ ■ • , AiAr))e 



e(^'Ar(/ii, • • • ,/iAr)|^'Ar(^i, • • • ,/UAr)), 
V .(^iv(M)IV'(xi,^l)l^7V+l({A})). .(^7V+l({A})|V;t(x2,t2)|^jv(M)). .... 

,i.tr A^N{{f^mNm))e.{^N+i{{xmN+i{{x}))e ■ ^ ^ 

aii|Aj-jv+l 



Using the equation ( |40| ) and the following relations : 

,(^jv(M)iv(xi,ti)i^^+i({A})), =, {^N+i{{\mHxi,ti)\^Nm))t, (42) 

,(M/;v({A})|M/,v({A}))e = (2L)^, (43) 
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we obtain 

X e^(*^-*^)S.- ,(*^+i({A})|V't(xi)|M/^(M)): ,(M/^+i({A})|Vt(x2)|*^(M))e. 
The translation invariance of time holds. 

,(^jv(M)|V'(3;i,il)V'^(x2,t2)|^Jv(M))e =e (^Jv(M) I V'C^^l , 0) (^2, ^2 " ^l) I ^-iV (M)). • (45) 

In the sequel we set the abberiviation t = t2 — ti. Remember a following simple fact. 

Lemma 4.1 For sequences {fji,--,j„}ji,--,j„el, {9ji,--,jn}ji,--,jn<^l^ ■ some index set), the follow- 
ing holds. 

J2 (Sym /)ji,-,in(Sym 5r)j,,...j„ = ^ (Sym g)j,,...j„. (46) 



Here we have used 



I 

(Sym f)j,,..,j„ = -J I] (47) 



The form factors have the determinant formulae in lemma 3.2 and 



N 

,(^^+i({A})|V^t(^)|^^({^}))^= J2 C7,(x|A^(^+i))n^.(^|A.o-),/i,), (48) 

o-e5jv+i i=i 



Therefore lemma 4.1, we obtain 



e 



^^^'^^-^^^^f^)'""^' E ••• E (e^*'--C:(-i|A^^-i)Q(x2|A^+0 + ^ 



X det (e**'^j/e(xi|Aj,/ifc)/e(a:2|Aj,/ij) - Qe**'^jC*(xi|Ao)Ie(x2|Aj,;Uj) (49) 

l<j,k<N \ 

X e**^^+iC7,(x2|A7v+i)/.(xi|A^+i,/ifc)) 

/ 0=0 

The j th line of the above matrix only depends on Xj not on A^, (A; ^ j), therefore we can insert the 
summations J2\i ' ' ' X^Ajv+i matrix. Now we arrive at the following. 

Proposition 4.2 The vacuum expectation values of two fields have the determinant formulas. 

,(^Af(/il, • • • ,/iAr)|V^(xi,0)-^t(x2,t)|^Af(//l, • • • , jUAf)), 
e(^'Ar(^l, • • • ,/XAr)|^'Ar(/Xi, • • • ,/UAr))e 
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ith 



2L 
1 f 1 



^ ^ ee^'''c,{xi\s)a{x2\s)+ ^ 



aa / i<j,k<N 



1 



— ae 



^ I ^ E e^*' C.(xi|s)Je(x2|s,/Xj) I I ^ I] e**' C7,(x2|s)J,(xi|s,^fc 



2L \ 2L 

Here we have used 

and functions Ce{x\s) and If:{x\s, fi) are defind in ([3^ and respectively. 



(50) 



a=0 



(51) 



The size of the above matrix depends on the state number A*", however, the element of the matrix 
does not depend on A^. By calculations, we obtan 

1 



2L 



ee**^'a(xi|s)Ce(x2|s) 



«efN 



2L 



E 



+ e— ^ gits ~is{x^+x2)^ (52) 
2^ «ef z 



— ^ e^*"'Ce(a;i|s)Je(x2|s,/i) 



(53) 



e 2 



and 



2 /tj-N ^its^-isxi 



vr 



E 

sefz 



s — fl 



sin(x2(s-/x)) -e^*^'-^/^^! ^ +e{ ^ ^ (-^) } 



e**'^ J,(xi|s,/i)J,(x2|s, A) 



(54) 



2 /vr 



2 /vr 



vr VL 



E 



1 + M(l + '^M,0 



\ — fl 



[e''^' sin(xi(A - ^)) + e^*^' sin(x2(A - //)) 



,S — A S — ;U 

Therefore we can take the thermodynamic limits. Let us set 

r(s, x) = its^ — ixs, 



sin((s - fi)xi) sin((s - A)x2) } + e—^ — { fi ^ (-^)} 

A + /i 



G{x) 



(55) 
(56) 



The following are our main results. 



Theorem 4.3 In the thermodynamic limit N,L ^ oo, such that = D, the ground state 
dynamical correlation functions have the Fredholm determinant representations. 



{ll}{xi,Q)^Hx2,t2))efi 



-iht 



2vr 



d 



G{xi - X2) + eG{xi + X2) + ^ ) det ( 1 - -V, - aW, 



vr 



0=0 



(57) 
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Here the integral operators and are defined by 

{vj){X)= rv,{x,fi)f{fi)dfi, {wj){X)= ^w,{x,^I)f{^^)d^l, (58) 





where the Fermi sphere q = ttD and the integral kernel are given by Neumann or Dirichlet sum : 

F,(A,^) = e-i^*("'+'^')(-^L(A,^) + e-^L(A,-/i)j, (59) 

I A - A + /i J 

W,{X,^i) = ee-^^*(^'+'''){P(A|xi,X2)+eP(-A|xi,X2)}{P(/i|x2,xi) + eP(-/x|x2,xi)}. (60) 
Here we have used 

L{X,n) = e**^%in(xi(A-/i)) +e^*'^%in(x2(A-^)) (61) 

2 ( \ 1 \ ■ 2 

H — / r e'*'' sin((s - ^)a;i) sin((s - A)x2)ds. 

vr 7-00 \s- II s - XJ 



( 1 


V. 




2 


POO 


vr 


J—oo 



P{X\xi,X2) = e"(^'^^)--/ -^e^^''^'hm{x2is - X))ds. (62) 

vr J-oo s — X 

We have suceeded to write the integral kernel by elementary functios : 

ds, (63) 

-oo s - y 

and trigonometric functions. Next we study the finite temperatute thermodynamics. By statistical 
mechanics arguments, at temperature T > 0, the thermodynamic equilibrium distribution of the 
spectral parameters is given by the Fermi weight 'i?(A) : 

1™ (t) \ \ =^(^^)' ^^^^^ ^(^) = WtY- (64) 

\LJ Xj+i - Xj 1 + exp ( 

Therefore the density is given by 

Dm = 

L vr 



N 1 

D{T) = — = - I ^{X)dX. (65) 



We arrive at the following theorem. 

Theorem 4.4 In the thermodynamic limit : N,L ^ oo,such that ^ = D(T), the finite temper- 
ature dynamical correlation functions have the Fredholm determinant representations. 



(V'(xi,O)'0"^(x2,t2)>e,r 



e 



-iht / Q 



G{xi - X2) + eG{xi + X2) + ^ det 1 V^^t - aWe,T 



(66) 



0=0 



27r V da J \ tt 

Here the temperature T > and the integral operators Ve^T md We.T are defined by 

(K,t/)(A)=/ v,{x,fiMf,)f{fi)dfi, {w,,Tf){X)= w,{x,^l)'^{^l)f{fl)d^l, (67) 



where the Fermi weight 'd^X) is given in (64)- 
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5 Degenerate Case 



The purpose of this section is to derive the Predholm minor determinant representations for finite- 
temperature fields correlation functions : 



{'4'{xi)'4'^X2))e,7 



(68) 



Especially our Fredholm minor determinant representations coincide with the one which has been 
obtained in O. When we take the limit t ^ 0, the following simplifications occur : 



-^x^ X 



G{x) ^ 0, L(A, /i) ^ sin(3;i(A — /i)) + sin(x2(A — /i)), P(A|xi, X2) ^ 
Therefore we obtain 

(V'(xi)^(x2)).,r = ^ [r^j det - -Ve,T - aW 



where the integral operators are given by 

^^.,t/) (A) 
and the integral kernels are given by 

v;,T(A,/i) 

+ e 



\ / -'0 



H^5/^(A,/x)/(^)d/i. 



^(A) 
1 



1 



{sin(xi(A - /i)) + sin(x2(A - ^))} 



_ A — /i 

{sin(a;i(A + /i)) + sin(x2(A + /u))} 



Pay attention to the Fourier transforms : 



/(A) 



1 



The following identity holds for functions /e(eA) = e/(eA). 

r dfiUfi)J^{X)^{fi) |^sin(x(A - ^,)) + e-^ sin(x(A + m)) 
Jo ^ lA-yU X + fi 

where 



e,T{C,ri)= / i?(i^){cos((C-?7)z^) +ecos((C + r/)zy)}(ii/. 
JO 



(69) 

(70) 
(71) 

(72) 

(73) 
(74) 
(75) 



(76) 



(77) 
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Therefore we arrive at 

de.(l-lf.,r)=det(l-^(«t/""')), 
where the integral operator 9^^^'^^^ is defined by 

^ ' ' Jyi 

Let us set 



(78) 



(79) 



det 1 



IT \ I 



Ae(C,r?) = 



det(l-f(C/^'^^0) 
Here we have used the following notation of the Predholm minor determinants : 

{-xr 



(80) 



det 1 - XKi 



m ■■■ Vr 

where we have used 



= E 



n=0 



6 



& Ai 

Vr Al 



Xr, 
Xr. 



,(81) 



6 



det {K{Cj,r)k)). 

l<j,k<m 



(82) 



m ■■■ Vm 

The integral operator Kj is defined by using the integral kernel K{X, jx) and the integral interval / : 

(83) 



(:^,/)(A) = J^K{X,fi)f{fi)dfi. 



Prom definition, the function Ag(^,77) satisfies the integral equation : 

2 p2 2 



2 r^^ 2 
A.(C,r?) - - / e,,T(^,e')Ae(e',??X = —9e,T{^,v). 

TT J -XI TT 



Let us take the Fourier transforms of this integral equation. 
1 



27rv/?(A) 



(84) 



(85) 



/ dfx%T{X, fi) ^ ^ / dC'e'f^^ A,(e', ry) = — 0?(A)(e^^'' + ee"*^"). 



Therefore we know 



A.(^,r?) = — TV 

TT 



TT 



(86) 



Now we arrive at 
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Corollary 5.1 The finite temperature field correlation functions have the Fredholm minor deter- 
minants representations. 



{^{x^)^Hx2))e,T = ^det ( 1 - 



X2 
Xi 



(87) 



where the integral operator 9^^^^'^^^ is defined in (7^) 



In much the same way as with finite temperature case, we obtain the Fredholm determinant repre- 
sentations for the ground state correlation functions. 

Corollary 5.2 The ground state correlation functions have the Fredholm minor determinant 
formulae. 



(V(xi)Vt(x2)).,o = Jdet ( 1 - -^(-1.-2) 
2 \ vr 



X2 
Xl 



Here the integral operator is defined by 

' Xl 

where 



^ ' J Xl 



Ke{i,v) = 7 + e TT • 90) 



Here the density D = ^. 



This result coincides with the one which has been obtained by the help of fermions Q. 
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